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Abstract 

A framework, which combines mathematical analysis, closure theory, and phenomeno- 
logical treatment, is developed to study the spectral transfer process and reduction 
of dimensionality in turbulent flows that are subject to rotation. First, we outline a 
mathematical procedure that is particularly appropriate for problems with two dis- 
parate time scales. The approach which is based on the Green’s method leads to the 
Poincare velocity variables and the Poincare transformation when applied to rotating 
turbulence. The effects of the rotation are now reflected in the modifications to the 
convolution of nonlinear term. The Poincare transformed equations are used to obtain 
a time-dependent analog of the Taylor-Proudman theorem valid in the asymptotic limit 
when the non-dimensional parameter y = Lit —>■ oo (SI is the rotation rate and t is 
the time). The ‘split’ of the energy transfer in both direct and inverse directions is 
established. Secondly, we apply the Eddy-Damped-Quasinormal-Markovian (EDQNM) 
closure to the Poincare transformed Euler/Navier-Stokes equations. This closure leads 
to expressions for the spectral energy transfer. In particular, an unique triple velocity 
decorrelation time is derived with an explicit dependence on the rotation rate. This 
provides an important input for applying the phenomenological treatment of Zhou. In 
order to characterize the relative strength of rotation, another non-dimensional number, 
a spectral Rossby number, which is defined as the ratio of rotation and turbulence time 
scales, is introduced. Finally, the energy spectrum and the spectral eddy viscosity are 
deduced. 
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1 Introduction 


Turbulent flows subject to solid body rotation have many important applications in aerospace 
engineering and geophysics. Although all turbulence models rely on assumptions on energy 
transfer, there is a great deal of difficulties in studying the time evolution of kinetic energy 
and dissipation of rotating turbulent flows. Although it has been long recognized that ro- 
tation can significantly modify properties of fluid flows, rotation does not even enter the 
transport equations for kinetic energy. This ‘indirect’ but significant influence by rotation 
is the physical root of the difficulties. A turbulence model which can account for rotation 
adequately has not yet been developed. 

Experimental studies 1-4 , direct numerical simulations (DNS) 5-8 , large-eddy simulations 
(LES) 5,9 , and closure approximation 10 have established that solid body rotation suppresses 
the nonlinear energy cascade from large to small scales. The effect of rotation is through 
phase scrambling. Zhou 11 recently noted the strong similarity between magnetohydrody- 
namic (MHD) turbulence and initially isotropic turbulence subject to rotation. For MHD 
turbulence, Kraichnan 12 pointed out that the propagation of the Alfvenic fluctuations dis- 
rupts the phase relation and thereby may be expected, on the average, to decrease energy 
transfer. Similarly, uniform rotation causes plane waves to propagate with phase speed 
2klk^/k. Here the rotation vector is considered to act along the vertical (|^y)- This evidence 
formed the basis for extending the phenomenology originally developed for MHD turbulence 
to rotating turbulence 11 . 

In this paper we focus our attention on developing a better understanding of the energy 
transfer process in homogeneous, isotropic, incompressible turbulence subjected to a uni- 
form background rotation. In particular, we address the problem of suppression of energy 
transfer and reduction of dimensionality due to rotation. We also deduce an energy spec- 
trum and spectral eddy viscosity that have an explicit dependence on rotation rate. Our 
analysis is based on the fact that this is a two time scale problem. We first present a general 
mathematical framework based on the Green’s method which is applicable when there are 
disparate time scales. The approach leads to the Poincare transformation and the Poincare 
velocity variables. The latter is accomplished with the time-dependent change of velocity 
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variables 


U(/,z) = L(Z)u = j G(nt,x,y)u(t,y)dy. (1.1) 

Here L(/) is the linear propagator describing propagation of inertial waves in rotating fluids 
(Poincare 13 ) and x,y) is the Green’s tensor corresponding to the gyroscopic operator. 

The effects of the rotation are now reflected in modifications of the convolution in the 
nonlinear term. 

The transformed equations are used to obtain a time-dependent analog of the Taylor- 
Proudman theorem valid in the asymptotic limit fi = 0/ — ► +oo. The subtlety of the un- 
steady situation lies in the fact that purely 2D turbulence is unstable to 3D perturbations 
and is, therefore, un-realizable. Bershadskii, et al. 14 argued that the quasi-2D turbulence 
regime arises as a result of a spontaneous breaking of reflexional symmetry, which in turn is 
a consequence of the instability of 2D turbulence to 3D wave motions. The experiments of 
Hophnger et al. 15 beautifully demonstrated that two-dimensional rotating fluid columns are 
capable of supporting wave motions. The major problem in obtaining a faithful description 
of quasi2D turbulence is to account for 3D effects. To this end, we deduce an analog of the 
Taylor- Proudman theorem for time-dependent flows. The operation of vertical averaging 
applied to the Poincare transformed equations is crucial in our analysis. The main conclu- 
sion is that in the above limit the vertically averaged velocity U satisfies the classical 2D-3C 
Euler/Navier-Stokes equations with initial data averaged in the direction. The difference 
V = U — U is, however, not small. The three-dimensional (^ 3 -dependent) vector held V 
is a passive vector transported by U. It satisfies equations with coefficients depending on 
U. The ‘split’ of the energy transfer in both direct and inverse directions is established: 
the vertically averaged velocity U transfers energy to large scales (inverse cascade) as 2D 
turbulence; and the passive vector V transfers energy directly from large to small scales. 

An EDQNM analysis 16,17 is then performed on the Poincare transformed rotating Navier- 
Stokes equations. Unlike the previous EDQNM analysis (EDQNM1 and EDQNM2) of Cam- 
bon and Jacquin 10 which was based on the original rotating Navier-Stokes equation, our 
analysis leads to a unique time scale for the decorrelation of triple velocity correlations. 
This provides an important input for applying the phenomenological treatment of Zhou 11 . 
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In this analysis, 73 , the time scale for decay of triple correlations, which is responsible for 
inducing turbulent spectral transfer, may depend on any relevant turbulence parameters. 
The usual assumptions, energy conservation by the nonlinear interaction and a local cas- 
cade, imply that e is independent of wavenumber k. Here e is the energy flux in the inertial 
range which also equals the dissipation rate. Local cascade implies that e is explicitly pro- 
portional to 73 and depends on the wavenumber and on the power of the omni-directional 
energy spectrum (Kraichnan 12 ). A simple dimensional analysis leads to 

e = A 2 T3(k)k 4 E 2 (k) (1-2) 

where A is a constant. The Kolmogorov 18 spectrum is recovered when T^ik) is chosen 
as a nonlinear time scale, T n i(k). Using the triple decorrelation time scale derived from 
EDQNM 16,17 , we deduce the energy spectra and spectral eddy viscosity at asymptotic limits. 
In order to characterize the relative strength of rotation, a spectral Rossby number, defined 
as the ratio of rotation and turbulence time scales, is introduced. For rotating turbulence 
in equilibrium, the Rossby number is the only relevant parameter controlling the effects 
of rotation on the flow. However, for the non-equilibrium situation, two non-dimensional 
parameters are required. 

Our analysis suggests the following picture of the energy transfer process in the limit 
of kit — >■ 00 and small spectral Rossby number (strong rotation). There are inverse energy 
transfers by U with —5/3 omni-directional energy spectrum. This is the same as two- 
dimensional turbulence. At the same time, there is also a direct energy cascade governed 
by the equation for the passive vector V at kit — ► 00 and 0 >> 1. The energy spectrum of 
the total velocity is c k~ 21 by Zhou 11 . 

This paper is organized as follows. In Section 2, a mathematical procedure based on 
Green’s method is described. The Poincare transformation which leads to the Navier- 
Stokes equations for the Poincare velocity variable is developed. The analog of the Taylor- 
Proudman theorem for time-dependent flows is derived in Section 3 by vertical averaging of 
the Poincare transformed equations. Analysis of the Poincare transformed energy transfer 
equations is presented in Section 4. The energy spectrum and spectral eddy viscosity are 
deduced in Section 5. The possibility of extension of the Green’s method to study the effects 
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of solid body rotation on flows in complex geometries is indicated in Section 6. 


2 Poincare transformed Euler/Navier-Stokes equations for 
rotating fluids 

We consider flows uniformly rotating with a constant angular velocity about the axis. 
The Navier-Stokes equations referred to rotating coordinates 

<9 

— U - z/V 2 U + U • VU + 20 JU = -Vo, (2.1) 

at 

V -U = 0. 


Here U = (Pi, P27 P3) is the velocity held, J is the rotation matrix, 2fi JU is Coriolis force 
and p is the modified pressure. We note that {JU}„= {e% x U}„= e„ 3 jUj where is the 
antisymmetric tensor of order 3, is not equal to zero only if ijl are different, and equal to 
1 if the permutation ijl is even and —1 if the permutation is odd. 

Eq. (2.1) can be written in the component form 

Ua,t T UjU at j T ^klc a 3 jU j — P,a (2.2) 

U a ,a = 0 . 

We introduce the Leray projection operator on the divergence free vector fields denoted by 
P. Applying P to Eq. (2.1), we obtain 

<9 

— U - z/V 2 U + P{U • VU} + 2DPJU = 0. (2.3) 

We note that P{U-VU}= -P{U x curl U}. The pressure is eliminated in favor of the tensor 
P(k) = {P a/ g(k)} (P a/ g(k) = S a [ 3 — k °‘ k 2 13 , k = |k| in Fourier space) using the incompressibility 
condition. Direct calculation shows that P(k)JU(k)= P(k)JP(k)U(k)= || e almkiU m (\i). 
Transforming to wavenumber space, the equation of motion (2.3) written in the component 
form reads 

o lx i 

— U a (k,t) + vk 2 U a {k,t) + 2kl-^e a i m kiU m (k) = -^P„/ m (k) ^ Ui(p,t)U m (q,t) (2.4) 

p+q=k 
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where 


Paf5",{ k) = kpP ai { k) + A; 7 P„ /3 ( k). (2.5) 

For fi = 0 we obtain equations which are the starting point of EDQNM analysis of non- 
rotating turbulence presented by Orszag 16 . 

Next we introduce the Poincare velocity u related to U as follows 

U = exp(-2fZPJZ)u = J G( — kU,x,y)u(y)dy. (2-6) 

Here G is the Green’s function of the gyroscopic operator (e.g. Ladyzhenskaya 19 ). The 
identity Eq. (2.6) can be easily inverted: u = exp(2fZPJZ)U = f G(flt, x, y)u(y)dy. We 
find that the Poincare velocity u is solenoidal: V • u = 0. The operator exp( — 2fZPJZ) 
is the linear propagator studied by Poincare 13 . The mathematical theory of the Poincare 
propagator has attracted a considerable amount of attention since the work of Sobolev 20 . 
In particular, a number of results on the structure of this operator in bounded domains are 
available (Ladyzhenskaya 19 ). They can be used to study the effects of background rotation 
on turbulence in complex geometries. In this paper we restrict our attention to homogeneous 
flows. Possible applications of the Green’s function method in complex geometries will be 
indicated in Section 6. 

Eq. (2.3) written in Poincare u-variables becomes 

d 

— u — z/V 2 u = exp(2fZPJZ){exp( — 2fZPJZ)u x cnr/(exp( — 2fZPJZ)u)}, (2-7) 

V • u = 0. 

The Poincare transformation (2.6) which relates velocity variables U and u can be 
written in the component form as follows 

U a ( k, Z) = cos(y2fH)n„(k, Z) - ^ sin(y2fZZ)e„ ;m /qn m (k, Z). (2.8) 

This transformation eliminates rotational term in (2.4) but alters the nonlinear term as 
shown below. 

It is convenient to introduce a 3 x 3 orthogonal matrix M = {M ai } defined on the linear 
subspace of vectors u (k) satisfying k • u(k)=0: 

M„ 7 (k, kit) = cos(^-2 klt)8 ai + j sm(^-2kU)e a i 1 ki. (2-9) 

A/ A/ A/ 
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The inverse of the 3x3 matrix M denoted by M 1 = {M ai } is easily computed as 


M( k, f It) = M (k, f It) = M( k, — f It). 


( 2 . 10 ) 


Then the Poincare transformation (2.8) is inverted as follows 

u a (k,t) = cos(y2fIf)U„(M) + t sm(^2klt)e a i m kiU m (k,t). 

Using the relation (2.8), we obtain 

c) k 

— U a (k,t) + vk 2 U a {k,t ) + 2kl-^e a i m kiU m (k,t) = 

cos(y2flt)^X(k, /) - t sin(y2f2t)e„ ;m A; ; ^y (k, t). 


( 2 . 11 ) 


( 2 . 12 ) 


Therefore, 

TO 1 / O 

cos(y2flt)XX(k,t) - -sin(y20/)e a ; m fc;-^(k,/) + vk 2 M am (k, -klt)u m (k, t) = 


-^Palm(k) Ul(p,t)U m (<l,t)(2.13) 

p+q=k 


Recall that the velocity vector U is expressed in terms of the Poincare velocity u via the 
Poincare transformation. 

Eq. (2.13) can be written in the equivalent form 

l -Pl/ m (k) ^ 

P2lm( k) 


M(k,-kU)( du ^^ + z/fc 2 u(k,f)) = - l - J2 Ui(p,t)U m (q,t) 


p+q=k 


\ P?,lm{ k) J 


(2.14) 


Inverting the 3x3 orthogonal matrix M( k, — kit) ( M x (k, —klt)=M( k, f It) by (2.10)), we 
obtain 


+ vk 2 u(k,t) = - l - Ui(p,t)U m (q,t)M(k,ttt) 

p+q=k 

We rewrite (2.15) in the component form 


1 Pllm( k) ^ 
P2lm( k) 

^ P?,lm{ k) ) 


(2.15) 
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cos(y20/)P a ; m (k) + sin(y2 ttt)e aia ^P a i m (k) 


—u a (k,t) + uk 2 u a (k,t ) 

Ui(p,t)U m (q,t). (2.16) 

p+q=k 


Substituting expressions for Ui and U m in terms of u via the Poincare transformation (2.8), 
we obtain the Poincare transformed Navier-Stokes equations written in terms of the Poincare 
velocity u as follows: 

d i 

—u a (k,t) + vk 2 u a (k,t)=-- ^2 F al 3 1 (k,p,q,2nt)u l 3(p,t)u 1 (q,t) (2.17) 


p+q=k 


where 


FaP'yik, p, q, 20/) = I] F 3 a p { k, p, q, 2 nt) = J2 / J ( k > P> q> (k:, p, q). (2.18) 


3 = i 


3 = 1 


Here 


QW k ’P>q) 


<5a/3 7 ( k 5 p 5 q) 

— Pa.(3m (k) ~€rnh'y Qh 

<5a/3 7 ( k 5 p 5 q) 

— Pal^(^P)~£lsf3Ps 

Qip-yO*-, P 5 q) 

— Palm(^-)~^lspPs~^mh-yQhi 

p q 

<5a/3 7 ( k 5 p 5 q) 

— ^a^cr ~j2'Po'/3'y(}P)j 

<5®/3 7 ( k ,p,q) 

hp /in <ih 

— ^a^cr j * a(3m \K-)€ m h'y i 

<5l/? 7 ( k 5 p 5 q) 

kt p n , Ps 

~ k " a F\ k ) e lsP •> 

<5a/3 7 ( k 5 p 5 q) 

k( -P Orl Ps 

— i -L<rlm\K-)€ls{3 ^mh'y 

K p 


Oh 

q 


(2.19) 


and 


f 3 (k, p, q, 20/) = ccc(k,p, q, 20/) = cos(-^20/) cos( — 20/) cos(— 20/), 


P 
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k* 


/ s (k, p, q, 20/) = sss(£;,p, q, 20/) = sin(^-2f2t) sin( — 20/) sin( — 20/). 


- 93 , 


A; 


P 


In the limit fi = 0/ — ► 0, the first term, j = 1, is the leading order term. The non-rotating 
limit can be recovered from the / 1 (k, p, q, 2 fit) term by setting 12 = 0 in Eq. (2.17): 


d i 

—u a (k,t) + vk 2 u a (k,t) = --P a i 3l (k) ^2 u p(P^ t ) u i( c l^ t )- 


p+q = k 


( 2 . 20 ) 


We note that the Coriolis term 2f2|fe„/ m /q?7 m (k) in Eq. (2.4) (after projection on the 
divergence free subspace) has the form of a scalar times the curl operator. Therefore, eigen- 
functions of this operator coincide with eigenfunctions of the curl operator. Decomposition 
of velocity fields along eigenfunctions of the curl operator is known as the “helical wave 
decomposition” and has attracted a considerable amount of attention since the work of 
Moffatt 21 , Moses 22 , Lesieur 17 and others. More recently, it was used in investigations of the 
energy transfer in both non-rotating and rotating turbulent flows 10,23 and for weakly com- 
pressible turbulence 24 . The main focus of these analyses is amplitude equations describing 
helical wave interactions. In the present paper we prefer to work directly with velocity 
variables without reverting to the helical wave amplitude equations. Our EDQNM develop- 
ment, presented in detail in the Appendix, is a natural extension to turbulence subjected to 
a uniform background rotation of Orszag’s 16 EDQNM analysis of non-rotating turbulence. 


3 Dynamics of vertically averaged velocity and the anal- 
ogy of the Taylor-Proudman theorem for time-dependent 
flows 

In this section we use the Poincare transformed equations to deduce an analog of the Taylor- 
Proudman theorem for time-dependent flows valid in the asymptotic limit fi = kit — ► +oo. 
The Taylor-Proudman theorem in its classical form states that steady inviscid flows in a 
rapidly rotating frame are two-dimensional, i.e. are independent of the coordinate along the 
axis of rotation of the fluid (Tritton 25 , Speziale 26 ). The importance of this statement stems 
from its use in calibration of turbulence models for rotating flows. According to Speziale 26 , 


8 



it is essential that turbulence models be developed whose properties in non-inertial frames 
of reference are consistent with the Taylor-Proudman theorem describing the structure of 
solutions of the Navier-Stokes equations in a rapidly rotating frame. 

The subtlety of the unsteady situation is, in the fact, that purely two-dimensional turbu- 
lence is unstable to three-dimensional perturbations and, therefore, unrealizable. The fact 
that two-dimensional rotating fluid columns are capable of supporting 3D (a^-dependent) 
wave motions is beautifully demonstrated in the experiments of Hophnger et a/. 15 . The 
major problem is to account for the collective contribution to the dynamics made by wave 
motions. To this end, we deduce an analog of the Taylor-Proudman theorem for time- 
dependent flows. 

The main idea is to obtain equations for the vertically averaged velocity using the for- 
mulation of the problem in terms of Poincare velocity variables. We find that the vertically 
averaged equations have a particularly simple form in the asymptotic limit f It — ► +oo. 
The main part of the equations coincides with the two-dimensional, three components (2D- 
3C) Euler/Navier-Stokes equations. The additional term involves terms varying only on 
a fast time scale. The structure of these equations allows us to obtain the analog of the 
Taylor-Proudman theorem for time-dependent flows. 

In this section we use the notation u a for the vertical (the direction of rotation) average 
of u a : u a = f u a dx 3 . Let k' = (Aq,&2,0) where k' is the 2D horizontal component of the 
wavevector k. Vertical averaging in Fourier space corresponds to restricting Fourier sums to 
wavevectors satisfying the condition = 0. Thus, n„(k) = M Q ,(k / ) is a 2D-3C velocity held. 
An important observation is that U=exp( — 2DPJt)u= u since M(k, 0/)|fc 3= o = Id (3 x 3 
identity matrix) which follows from (2.9). Applying the operation of vertical averaging 
to Eq. (2.16), we obtain 


d i 

—u a (k',t) + isk 2 u a (k',t) = --P a i m (k) Ui(p,t)U m (q,t) = 

p+q=k 


(3.1) 


£ 

p+q=k,fc 3 =p 3 +g 3 =0 


Recall that the velocity vector U is related to the Poincare velocity vector u by Eq. (2.8). 
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Then Eq. (3.1) together with k% = 0 implies 


dt 

where 


d i 

—u a (k , ,t) + vk 2 u a (k , ,t) = -- Y Th /37 (k / ,p,q,2fit)M /3 (p,t)M 7 (q,t) (3.2) 


p'+q'=k',p3+g 3 =0 


4 4 

£*/? 7 (k',P,q,2ftZ) = = ^/ J ( k ',P,q, 2 ^)Qa^ 7 ( k, >P>q)- ( 3 - 3 ) 

3 = 1 4 = 1 

We stress that k% = 0 implies that we have only four terms in the summation. Since 
k 3 = p 3 + q 3 = 0, we obtain 


/ 1 (k, p, q, 20/) = cos(— 20Z) cos(— 20Z) = cos(— 20Z) cos( — 20Z), 

p q p q 

/ 2 (k, p, q, 20/) = cos(— 20/) sin(— 20/) = — cos(— 20/) sin(— 20/), 

p q p q 

/ 3 (k, p, q, 20/) = sin(— 20/) cos(— 20/) = sin(— 20/) cos(— 20/), 

p q p q 

/ 4 (k, p, q, 20/) = sin(— 20/) sin(— 20/) = — sin(— 20/) sin(— 20/). (3.4) 

In the hmit p = kit >> 1, Eq. (3.2) has a separation of time scales and can be further 
simplified. Since f 2 and f 3 are expressed as sums of sin functions only, these terms make 
no contribution to the dynamics on slow time scale (no resonances). On the other hand, the 
functions / 1 and / 4 are expressed as sums of cos functions. Then resonances contributing 
to the slow time scale dynamics are possible only in the case p = q. Combining these terms, 
we obtain 


/ X ( k', p, q, 2m)Q^ 7 (k', p, q)up(p, Z)u 7 (q, t) + 

/ 4 (k', p, q, 20Z)Q 4 /37 (k / , p, q)n /3 (p, Z)ti 7 (q, t) = 

= cos 2 (~20Z)P„ / 3 7 (k)n / 3(p, Z)n 7 (q, t) - sin 2 (^20Z)P„ /37 (k)n /3 (p, Z)w 7 (q, t) = 

7 ) o 

= cos(^4klt)P a p 1 (k)up(p,t)u 1 (q,t). (3.5) 

In the derivation of Eq. (3.5) we have used the fact that Q^p (k, p, q) = P a/ g 7 (k) if = 0, 
p = q. The above calculation shows that 

- - Y F arjl {k\ p, q, 2klt)u fj {p , Z)n 7 (q, t) = 

p'+q'=k',p3+g 3 =0 


Y P«/ 3 7 (k / )n / 3 (p / ,Z)n 7 (q / ,Z) + P(0Z, u, u) (3.6) 

p'+q'=k>3=g 3 =0 
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where R(Tlt, u, u) denotes terms evolving on a fast time scale only (without resonances 
contributing to the dynamics on the slow time scale). Finally, Eq. (3.2) becomes 

d i 

—u a (k',t) + nk 2 u a (k',t) = -- P„ / 3 7 (k / )n / 3 (p / ,t)n 7 (q / ,t) + P(0/, u, u)(3.7) 

p'+q'=k>3=g 3 =0 

where _R(0/, u, u) is a rapidly oscillating term 

1 /* T COT! st 

|- jf R(Slt, u ,u)dt\< — ^0 (3.8) 

for flT >> 1. Eq. (3.7) written in physical space has the form 

d 

— u - z/V 2 u + u • Vu = — Vp + R(Vlt, u, u), (3-9) 

V • u = 0. 

Thus, up to the rapidly oscillating terms f?(fP,u,u), the vertically averaged velocity 
u = U satisfies the 2D-3C Navier-Stokes equations (3.9) which are averaged in x% initial 
data. Therefore, we obtain 

The analog of the Taylor-Proudman theorem for time-dependent flows: 

The vertically averaged component U of the exact solution U of 3D-3C Euler /Navier-Stokes 
equations with initial data Uo is close to the exact solution of 2D-3C Euler /Navier-Stokes 
equations U with initial data Uo averaged in x 3. The difference V = U — U between the 
full 3D-3C velocity U and 2D-3C velocity U corresponding to the vertically averaged initial 
data is, however, not small. 

A mathematically rigorous formulation and proof of the analog of the Taylor-Proudman 
theorem for time-dependent flows (including error estimates) are presented in Babin, Ma- 
halov and Nicolaenko 27 . 

We now present an alternative development of the time-dependent analog of the Taylor- 
Proudman theorem based on physical arguments. Both computations 5-9 and experiments 1-4 
have noted an increase in integral lengthscales along the rotation axis relative to those in 
non-rotating turbulence. An increase in the integral lengthscales has been thought to be 
a prelude to a Taylor-Proudman reorganization to two-dimensional turbulence. It is rea- 
sonable to conjecture that at zero Rossby number the only active triads in the asymptotic 
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state are those which are invariant under these scalings. This is called the selection principle 
for triad interactions in the asymptotic state. Clearly, 2D-3C interactions corresponding to 
p 3 = q 3 = 0 are invariant under scalings of p 3 , q 3 in a trivial way. On the other hand, in the 
asymptotic state of rotating homogeneous turbulence interactions are restricted to pairs of 
wavevectors satisfying the resonance conditions 28,23 . 


± P3 ± 03 = h 
p q k 


(3.10) 


One can easily find solutions of the resonance equation (3.10) obeying the selection 
principle. The solutions are {^3 = 0 ,k = p}, {p 3 = 0 ,k = q} and {^3 = 0 ,p = q}. 
As shown above, the interactions {^3 = 0,p = q} make no contribution to the dynamics 
on slow time scale. Two other resonance conditions can be written in equivalent form 
{g 3 = 0, k\ + k 3 = Pi + P 2 } and {P 3 = 0, k\ + k 3 = q\ + q 3 }- Clearly, they are invariant 
under the process of an increase in integral lengthscales along the rotation axis relative to 
those in the horizontal plane ( k 3 — ► fk 3 ,p 3 — ► fp 3 ,q 3 — ► ^3). Thus the only active triads 
in the asymptotic state are {^3 = 0, k\ + k 3 = p\ + p 3 } and {p 3 = 0, k\ + k 3 = q\ + q 3 }- 
The condition p 3 = 0 simply means that the first entry in the bi-linear form is replaced 
by its averaged (vertical) value. Similarly, the condition g 3 = 0 implies that the second 
entry is replaced by its averaged value. Therefore, the scale-invariant equations valid in the 
asymptotic limit p = kit — ► +00 have the form 


<9u 

~di 


B(u, u) + B(u, u). V • u = 0. 


(3.11) 


Here u denotes vertical averaging of velocity. Applying the operation of vertical averaging 
to Eq. (3.11) and using the fact that B(u, u) + B(u, u)= B(u, u) + B(u, u)= — u • Vu, we 
conclude that u = U satisfies the classical 2D-3C Euler equations. 

The analog of the Taylor- Proudman theorem for time dependent flows provides a clear 
picture on the reduction of dimensionality at the limit kit — ► 00. The flow is three- 
dimensional when the R term in Eq. (3.9) is significant (no ‘split’ of velocity into U and 
V) while the flow is quasi-two dimensional when R is negligible (‘split’ of velocity into U 
and V). The error estimate, Eq. (3.8), can be used to obtain the value of critical p* when 
such a transition takes place. 
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We are now in a position to present a physical picture when the time-dependent flow 
satisfies the analog of the Taylor- Proudman theorem for time dependent flows. There are 
many two-dimensional layers. If vertically averaged, then the dynamic equations for U 
are governed by the 2D Navier-Stokes equations with fast-time scale corrections (Eq. 3.7). 
In this case, the direction of the energy transfer for U is from direct to inverse cascade 
within each pancake. These layers, however, are not isolated among themselves. For the 
Euler equations, the ^ 3 -dependent component V is not small even at zero Rossby number 
as shown in Babin, Mahalov and Nicolaenko 27 . The evolution equation for V satisfies a 
‘passive vector’ equation that is driven by U. This equation couples these two-dimensional 
layers and is responsible for the (suppressed) direct energy transfer from large to small 
scales. The energy spectrum in this case will be discussed in Section 5. 

4 The Poincare transformed energy transfer equation 

In this section the EDQNM closure is applied to the Poincare transformed Navier-Stokes 
equations written in terms of the Poincare velocity u. Following the standard EDQNM 
closure procedure, the triple Poincare velocity correlation is expressed in terms of second 
order moments. The details are presented in the Appendix. The novelty of our approach 
is that the EDQNM closure is applied to the Poincare transformed equations which shows 
the need for a heuristic approximation for the effects of rotation in the third order moment 
equations. Lesieur 17 points out that choice of eddy damping is more difficult in problems 
where waves interact with turbulence. The previous EDQNM analysis of rotating turbulence 
requires three evolution equations for second order moments [Cambon and Jacquin 10 ]. The 
Coriolis effect, fi, only appears in one of them. As a result, there are two possibilities on 
the memory integral. The EDQNM1 version of Cambon and Jacquin 10 does not include 0 
in the ‘damping rate’. On the other hand, 0 is introduced explicitly into in their EDQNM2 
version. They found that these two choices will lead to different results and the latter 
is needed for the anisotropization mechanisms and to predict a departure from isotropy. 
Our EDQNM closure is based on the Poincare transformed Navier-Stokes equations. This 
Poincare transformed equation has the same spectral structure as that of the original one 
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(with 0 introduced explicitly into the generalized projection operator) and has a unique 
memory integral. 

The EDQNM closure of the Poincare transformed equation, presented in detail in the 
Appendix, provides a framework to study the observed suppression of energy transfer when 
turbulence is subject to rotation. We assume isotropy of the energy spectrum when exam- 
ining the energy transfer process. Experiments 1-4 and DNS 5-9 of rotating homogeneous 
turbulence established that although the rotation increased the anisotropy of the length- 
scales, the three principal components of the Reynolds stress tensor remain nearly equal 
provided that the initial state is isotropic. Thus, our assumption is consistent with these 
findings. 

The isotropy assumption leads to 

^(k, t) = 7>U(k, t)P al3 ( k), (4.1) 

where U(k , t) = and E(k, t) is the energy spectrum. From the Appendix, the equation 

for U(k , t) reads 

+ 2 vk 2 )U{k,t) = ^(a'(A;,p,g,ff,t)E(p,t)E(g,t)+ (4.2) 

V 

b'(k,p, q, 0, t)U(k, t)U(q, t ) + c\k,p, q, 0, t)U(k, t)U(p, /)), 


where k = p + q is understood. The geometrically determined coefficients a'(k,p,q,Cl,t), 
b'(k,p,q,Cl,t) and c'(k,p,q,Cl,t) are given in the Appendix. We note that they contain / 
and the memoral integral 0 (see (A. 5)). In order to evaluate these integrals, we expand / 
in a sum of sin and cos using the identities of a kind 

cos (j) cos ip cos x = — (cos(<^> + "0 + x) + cos(<j) + ip — x) + cos(<p — ip + x) + cos(<p — ip — x)), 

Odd terms in t can be expanded in sums of sin functions. For example, 

sin cp cos ip cos x = — (sin(<^> + V* + x) + sin(<^> + "0 — x) + sin(<^> — "0 + x) + sin(<^> — ip — y)). 


These identities and the following identities 


e ax cos (bx)dx = 2 _[_ 52 ( a cos (^ x ) + & sin(6a;)), 

f e ax 

j e ax sm(bx)dx = -^——j^(asm(bx) — bcos(bx)) 
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are now used to obtain analytical expressions for the coefficients 0 (see (A. 5)). 

Let 6 = , ib = 2a. and y = US-. Then we obtain 

^ |p+q|’ T p A q 

n -J 

/ e ~\y( k +'P +9 )+Pkpq](t-r) cos ((()2ffr) cos(-02fir) cos(x2ffr)dr = -(0 ++ (A;,p, g, 0.0/) + 
do 4 

0_l (A;,p, g, 0, fit) + 6 g, 0, fit) + 0__(A;,p, g, 0, fit)), 


/ e _ M fc +p +9 )+w= P ?](*- T ) C os(<^2flr) sin(-02fir) sin(x2fir)dr = -(6 . | (A;,p, g, 0.0/) + 

do 4 

0 |_(£;, p, g, 0, fit) — 0_|_ + (A;,p, g, fi, fit) — 0__(A;,p, g, fi, fit)), 


/» ■£ -j 

/ e -[i / (fc 2 +p 2 +9 2 )+A‘fc Pq ](t-^) s in(())2fir) cos(rd2fir) sin(x2fir)(tr = -(0_| (A;,p, g, fi.fit) + 

do 4 

0__ (k,p, g, fi, fit) — 0_|__|_(£;,p, g, fi, fit) — 0 |_(A;,p, g, fi, fit)) 


/ e -[i / (fc 2 +p 2 +9 2 )+A‘fc Pq ](t-^) s in(())2fir) sin(rd2fir) cos(y2f It)c1t = -(0 (A;,p, g, fi.fit) + 

do 4 

0 |_(&,p, g, fi, fit) — 0_|__|_(A;,p, g, fi, fit) — 0_| (k,p, g, fi, fit)). 

The EDQNM relaxation time scales 0±±(&,p, g, 0, Ot) and 0±±(&,p, g, 0, Ot) are given by 




6±±(k,p,q,tt,ttt) = - 

(<f> ± ip ± x)2fi sin((<^> ± ip ± x)2f it) + [//(A; 2 + p 2 + q 2 ) + cos((<f> ± ip ± x)2f it) 


// 2 (A; 2 + p 2 + g 2 ) + Pfcpg] 2 + ± ip ± x) 2 fi 2 

[//(A; 2 + p 2 + g 2 ) + 


- exp [-[//(A; 2 + p 2 + g 2 ) + Pfc pg ]t] 


[// 2 (A; 2 + p 2 + g 2 ) + p^] 2 + (<£ ± -ip ± x) 2 fi 2 


}(4.4) 


and 




6±±(k,p,q,tt,ttt) = - 

(<f>± / ip ± x)2f i cos((<f> ± ip ± x)2f it) + [//(A; 2 + p 2 + g 2 ) + sin((^ ± ip ± x)2f it) 


[// 2 (A; 2 + p 2 + g 2 ) + p^] 2 + (<£ ± i> ± x) 2 fi 2 
■2,2, 2x , „ 1,1 (0±^±x)2fi 


- exp [-[//(A; 2 + p 2 + g 2 ) + n kpq \t\ 


[v 2 (k 2 + p 2 + g 2 ) + p^] 2 + (<£ ± ip ± x) 2 fi 2 


} (4-5) 
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Since the eddy damping is much greater than the molecular viscosity, the EDQNM time 


scale for the triple velocity decorrelation can be simplified to 


/) n„ „ „ O nn ~ , 1 r ± ^ ± x)2D sin((^ ± kp ± x) 2 Dt) + p kpq cos((^ ± ip ± x) 2 Dt) 


/ 4 P3 + (0 ± t /’ ± x ) 2 ^ 2 

- exp 


(^kpq 


Vkpq + 


X4-6) 


and 


h n„ „ „ O nn 1 r ± x)2D cos((^ ± ip ± x)2Dt) + Mfcpg sin(0 ± ip ± x)2Dt) 

V±± («;, P, (£, St, ‘*Zj ~ \ 2 . / I I I I \2Ci2 

4 Mfcpg + ± k ± X) 2 ^ 2 

r ( 0 ±^±x) 2 ft V/l ^ 

exp[ Pl " i /‘1„ + (<* ± ¥’ ± 0 2 ff-X 4 7) 

They characterize the relaxation of < u(k)u(p)u(q) > toward a quasi-equilibrium by 
nonlinear transfer and molecular viscosity. Rather than solving the EDQNM evolution 
equations, we will only use information about the triple time scale of velocity decorrelation 
obtained from this theory as an input for a phenomenological treatment. In the next section 
we will use this information to deduce the energy spectrum and the spectral eddy viscosity 
in various asymptotic limits. 


5 Phenomenological analysis at asymptotic limits 

In order to infer the form of the inertial-range spectrum, it is necessary to estimate the 
magnitude for the triple correlations. In general, 73, the time scale for decay of triple 
correlations which is responsible for inducing turbulent spectral transfer, may depend on any 
relevant turbulence parameters. Because energy is conserved by the nonlinear interaction 
and a local cascade has been assumed, e is independent of k. Local cascade 29-31 also implies 
that e is explicitly proportional to 73 and depends on the wavenumber and on the power of 
the omni-directional energy spectrum 12 . A simple dimensional analysis leads to 

e = A 2 T 3 (k)k 4 E 2 (k ) ( 5 . 1 ) 

where A is a constant 32 . If the energy- containing range excitation were absent, it would 
be expected that the nonlinear interaction would build up substantial triple correlations in 
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the local dynamical time T n i(k). In steady state isotropic turbulence, the only time scale 
available for r^{k) is the nonlinear time scale, T n i(k ); therefore the Kolmogorov 18 spectrum 
is recovered. The equation which determines the spectral eddy viscosity, vj(k), is given 
by 33 

e = v T {k)S 2 {k) = v T {k)E{k)k 3 . (5.2) 

In our spectral representation, it is convenient to use the spectral Rossby number 
Ro s (k) = /ifc/O, which is the ratio of turbulence (nonlinear eddy turn-over) time scale 
and rotation time scale. For a given inertial Rossby number, we can define the cases of 
strong rotation when Ro s <C 1 and the weak rotation when Ro s >> 1. 

We now consider the number of the non-dimensional parameters needed. For the turbu- 
lence in equilibrium the Rossby number is the only relevant parameter controlling the effects 
of rotation on the flow. However, for the non-equilibrium situation, new non-dimensional 
parameters fi = kit and a = are required. They correspond to non-equilibrium effects 
at long and short time scales, respectively. 

We have shown that the parameter n = kit is introduced into the time scale of the 
triple velocity product as the arguments of sin and cos functions. These terms are clearly 
responsible for the scrambling of the energy transfer and do not go to zero as t — ► oo. As 
a result, the values of the time scale of the triple velocity product will oscillate; and this 
oscillation can be removed only by a long time average. For a turbulent flow with a given 
initial Rossby number (therefore, the rotation rate 1/0), the non-dimensional parameter fj, 
will evolve in three phases: (1) fj, — ► 0; (2) fj, bounded; (3) fj, — ► oo. In addition, a = [i^t is 
introduced in the second term in both 0± and 0± as a damping factor. This damping term 
is significant initially but it will go to zero as / — ► oo. 

It is clear that the general treatment of this problem is rather complicated. When 
a resonance condition, say cj) + ip + x = 0, is satisfied, the corresponding decorrelation 
time 0 ++ is simply the nonlinear time scale T n \ = l/[ k 3 E(k)Y^ 2 with a damping factor 
[1 — exp [— for short time scale. 

The non-resonance condition needs to be examined with care. Asymptotic limits of 
strong and weak rotation offer considerable simplification. These limits correspond to small 
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(spectral or initial) and large Rossby numbers, respectively. 
For the small Rossby number case, we have 


Q±±(k,p,q,Sl,Slt) 


sin((<)> ± ip ± y)20/) 

(0±^±x)fi 


- exp [-Pkpgt] 


(^kpq 


(<f> ± p ± x) 2 n 2 


( 5 . 3 ) 


and 


0±±(k,p,q,n,tlt) 


cos((<?5>± ± y)20f) n 1 

(< p±k±x)n expL f * kpqi (<t>±i>± x )n 


Alternatively, for large Rossby number, one can show that 


( 5 . 4 ) 


a / / o on., C0s ((<(>±l/’±X)2^) r 4l 1 

v±±(k,P, q, SI, Sit) ~ exp [-fikpqt \ 


ftkpq 


ftkpq 


( 5 . 5 ) 


and 


a n O _ sin((<)> ± r/> ± x)2Dt) (f> ± if ± x)2D 

0±±(«bP> q, SI, Sit) & exp [-fi k pqt\ 2 • 

^kpq P k pq 

5.1 Disruption of homogeneous isotropic turbulence by rotation 


( 5 . 6 ) 


We now focus on the short time behavior of the decorrelation time scale for the triple 
velocity products. 

There are two different situations. First, the time under consideration is so short (t — ► 0) 
that both 9±± — ► 0 and 9±± — ► 0 by definition in Eqs. (4.6)-(4.7). When the rotation is 
also very strong, it is called the Rapid Distortion Theory (RDT) limit 34-36 . To make 
contact with RDT, we rewrite p = t* /t^, where t* = e t/I( 0 and = e/(K 0 £l). This 
procedure nondimensionalizes the time scales using the initial turbulent kinetic energy (Ko) 
and dissipation rate. Note that the inertial Rossby number of the flow is defined as Ro = 
0.5e/(ifo^)- RDT theory assumes that the nonlinear effects do not have time to build up 
under the conditions of very short time and strong system rotation. Specifically, the RDT 
limit is valid when the initial Rossby number Ro <C 1 and the time normalized by initial 
turbulent time f < 1. 
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We now consider the second case where the nonlinear effects become important but time 
t* is still short. Again, two distinct physical situations are included in this limit: (a) weak 
rotation; (b) strong rotation (we call it ‘disruption’ as an extension of the RDT limit). In 
the ‘disruption’ limit the rotation is very strong. In a fixed time interval, II is very large. 
However, the time t is very short so that Clt is bounded. 

For the weak rotation case (large Ro s ), we have 

9 ss 0 ~ [1 - x — — . (5.7) 

f^kpq 

Similarly, for the strong rotation case, one has 

0 ~ [l-e-^*] x sm @ nt \ (5. 8 ) 

and 

0 ~ [l - e -w=pA] x cos ^ nt ^ . (5.9) 

These relations form a basis for constructing implicit equations for the energy spectra 
and spectral eddy viscosity in non-equilibrium turbulent flows. Such flows are deferred to 
a future study. 

5.2 Long-time dynamics of uniformly rotating fluids 

We now consider the limit of long-time dynamics. Here, the long time limit means that 
the time t is large enough so that the damping terms in 9 and 9 can be neglected. Again, 
we must consider two physically distinct situations: strong and weak rotation. In the weak 
rotation case we find that the decay rate of triple velocity correlations is given by the 
turbulence nonlinear time scale T n {. Therefore, at leading order, the energy spectrum and 
spectral viscosity are the same as that of the non-rotating case. The rotation effects can be 
introduced trivially here by keeping the next order correction for rotation in 9 and 9. We 
omit this analysis here. 

In the the strong rotation case, the dependence of 9±± and 9±± on H should be consid- 
ered with care. The time scale for Ts(k), the decorrelation of the triple velocity product, is 
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again the controlling parameter influencing the energy transfer process. When a resonance 
condition, say (f> + ip + x = 0, is satisfied, the corresponding decorrelation time 0 ++ is simply 
the nonlinear time scale T n \ = 1 /[k 9 E(k)Y^ 2 ■ In the mean time, the other contributions to 
this time scale at this given condition, 0 . | , 0 | are of order 

0 ~ 0+_ ~ 0_+ ~ 0__ ~ sin (^ flf ) . (5.10) 


Similarly, we find that 




cos(20/) 

ft 


(5.11) 


In order to remove the phase from Eqs (5.10)-(5.11), we introduce the long time averaged 
time scale which determines the decorrelation of triple velocity products at a small spectral 
Rossby number ( ft/ — ► oo): 


r 3 = 0 Q = yj 91 + §1 ~ I/O. 


(5.12) 


Therefore, in the regime of high Reynolds numbers and low Rossby numbers, rotating 
turbulence is characterized by two disparate time-scales: a short time scale associated 
with the rotation frequency tq = I/O and a nonlinear time scale. Our EDQNM time 
scale discussed above provides an additional support for the choice of tq = I/O in the 
phenomenological analysis of Zhou 11 . This is the key assumption in his application of the 
extended phenomenology of Kraichnan 12 and Mattheaus and Zhou 32 to isotropic turbulence 
subject to rotation. This is also the basic assumption in the large-eddy simulation (LES) 
of strong rotating turbulence of Squire et al. 9 where they argued that the time scale of the 
triple velocity correlation is proportional to tq = I/O. Squires et al. 9 assumed that the 
correlation time of the nonlinear triadic interactions is directly proportional to the short 
time scale tq. The predictions of the asymptotic decay for the kinetic energy in rotating 
turbulence, based on this crucial assumption on the time scale, have been confirmed within 
a few percent by their LES runs for times of order O(10 3 ) of initial turbulence. We find 
that a direct application of 73 = tq resulted the energy spectrum for turbulence subject to 
strong rotation: 

E(k) = C Q (kle)^ 2 k- 2 . (5.13) 
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In general, Zhou 11 has pointed out that there is a need to incorporate both rotating and 
eddy-turn-over time scales. We make use of the viewpoint that the lifetime of triple corre- 
lations in rotating turbulence might be more accurately treated by taking into account the 
possibility that these correlations decay because of the influence of both wave propagation 
and nonlinear triadic interactions 32 . The simple choice 


Tzik) T n i(k) + m(k) 

satisfies the appropriate limiting cases: T^k) — ► r n i without rotation and T^k) — ► tq with 
strong rotation. The generalized inertial range energy spectrum is 


E(k) = Z 2 A~^ 3 ^ 3 k~^ 3 , 


(5.15) 


where Z is given by 11,32 


[VY+ \]-Y + 2/h 2 + 4Z 0 ) 


z= 1 -{Vy + 


(5.16) 


where 


Y = 


1 1 ,4Z(k, 

\j2 + Vi + ( — ,3 + 


\] I “ Vi + 


. 4Zn 


(5.17) 


The parameters A = (Matthaeus and Zhou 32 ), k q = (fl 3 /e) 1 / 2 and Zq = 

A 2 / 3 ^ fc2 ^i /3 = [“ f 2 "] 2/,S - The strong rotation limit then leads to C'q = 1/A = 1.22 — 1.87 
for the typical range of Kolmogorov constant. These equations reduce to the non-rotating 
Kolmogorov “-5/3” spectrum when Zq — ► 0 (so that Z — ► 1), and to our rotating modified 
“-2” spectrum when Zq — ► 00 (so that Z — ► ^o^ 4 ). Fig. 1 illustrates spectra obtained this 
way with e = 1 and A = 0.7 (we have chosen Ck = 1-5), varying the reference rotation 
rate over values kl = 0, 10, and 100. The zero rotation case is a pure Kolmogorov spectrum 
while the kl = 100 case is very nearly k~ 2 spectrum. For intermediate rotation rates the 
spectrum varies smoothly between these two limiting forms, according to the increase of 
the controlling parameter Zq with increasing ratio kn/k. 

The rotation dependent eddy viscosity can be estimated 11 as 

E(k) v m 1 


v T (k) = M^r 1 ] 1 ' 2 


1 + "M -1 / 2 ’ 


(5.18) 
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where vq is a constant and y = k 3 E(k)/kl 2 is defined following Zeman 33 . The inertial range 
wavenumber k can be related 37-38 to the turbulent kinetic energy K and dissipation rate e 


* = fw (3C '‘/ 2 > 3/2 ' 


For the inertial range spectrum, Eq. (5.18) can be rewritten in physical space: 


vt{x) = v' 0 


,K 2 


1 


(5.19) 


e 1 + 0.36 KSl/e’ 

The eddy viscosity above is only appropriate for homogeneous turbulence without mean 
velocity gradients. For turbulent flows with mean velocity gradients in a rotating frame, 
the effect of rotation appears in conjunction with the mean vorticity 39 and Eq. (5.19) can be 
extended in a straightforward way by replacing by ( WijWij ) lj/2 (Gatski and Speziale 40 ). 
Here W \j = Uij + ae m jikl m is the extended absolute vorticity tensor defined by Gatski and 
Speziale 40 (a is a constant). The constant v' Q can then be evaluated in homogeneous shear 
flow where Wij reduces to the uniform shear rate. Such flows are deferred to a future study. 
We note that the scaling of the eddy viscosity with rotation rate is inconsistent with RDT. 

The spectral time scale is an important measurement 11 . The Kolmogorov hypothesis 
implies that the energy- containing range excitation does not affect energy transfer within 
the inertial range. Therefore, the average rate of energy dissipation is identified with the 
rate of spectral energy transfer and the rate of energy input. The nonlinear (or eddy 
turnover) time-scale, T n i(k) = [k 3 E(k)\~ 1 l 2 ^ is then equivalent to the spectral transfer time, 
t s . Zhou 11 has shown that 

1 [ T nl(k )] 2 


or more generally, 


T.s(k) = 


T s (k) = 


A 2 Tn(k) 

1 T nl( k ) 

A 2 r 3 (k ) ’ 


(5.20) 


(5.21) 


Therefore, the time for the spectral energy transfer is increased to a value greater than T n i; 
thus nonlinear energy transfer is suppressed by rotation. 

We stress that the long-time averaging procedure here is the same as that introduced 
for the time-dependent analog of the Taylor- Proudman theorem. This averaging procedure, 
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as defined in (3.8), removes the oscillations and ‘splits’ the velocity fields into two parts. 
Our analysis suggests that the energy transfer process in the limit of Of — ► oo and small 
spectral Rossby number (strong rotation) is as follows. There are inverse energy transfer by 
U=u with —5/3 omni-directional energy spectrum. This is the same as the two dimensional 
turbulence. In the meantime, there is also a direct energy cascade governed by the equation 
for the passive vector V introduced in Sec. 3 at Of — ► oo and 0 >> Hk- The energy 
spectrum of the full velocity U is given by Eqs. (5.15)-(5.17). 


6 Concluding Remarks 


In this paper we presented a general framework to study turbulence subjected to a uni- 
form background rotation based on Green’s method. In the analysis presented above we 
assumed spatial homogeneity of turbulent hows. Although the application of this method 
in complex geometries alters the integral convolution kernels in the Poincare transformed 
Navier- Stokes equations, the mathematical procedure is still valid provided that explicit 
asymptotic expressions for the Green’s tensor can be obtained. For example, in the case of 
a turbulent how subjected to an uniform background rotation occupying a region D, the 
Green’s tensor corresponding to the gyroscopic operator is found from the linear problem 

d , 

— U - z/V 2 U + 2fiJU = -Vo, 

at 

V-U = 0. (6.1) 


Then the Poincare transformation which relates the fluid velocity U with the Poincare 
velocity variable u is given by 

XJ(t,x) = Gu = f G(Clt,vt,x,y)u(y)dy. (6-2) 

Jd 

As a result, the Navier- Stokes equations written in terms of the Poincare velocity variables 
take the form 


d 

— u = B(0/,z4, u, u), 


(6.3) 


where 


B(0/, vt, u, u) = G 1 P{G'u x curlG u}. 


(6.4) 
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When explicit expression for the gyroscopic operator G can be obtained, the relation be- 
tween the nonlinear (turbulent), rotation, and viscous time scales in Eq. (6.4) are trans- 
parent. In general, explicit expressions are not available. However, for certain domains of 
interest (e.g. turbulent flow between rotating parallel plates) one can obtain explicit expres- 
sions for the gyroscopic operator valid in the asymptotic limits p = fit — ► 0, p = fit — ► +oo. 
Then one can apply the analysis presented in this paper. 
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A Appendix: EDQNM closure for the Poincare transformed 
equations 

We consider the most general term in (2.17). It has the form 

^2_l 2 , 

~Qf(k ,t) + vk 2 u a (k,t) = --'^2f{p,k,M)Q afjl {k,p)uf j {'p,t)u 1 {k - p ,/). 

V 

We define 

S af3 (k,t) =< u a {kp)u fj {-kp) >, 

TafoiK P, t) =< M„(k, t)u fj { p, t)n 7 (-k - p, /) >, 

U af 3 lS ( k, p, q, t) =< u a ( k, t)u fj { p, t)n 7 (q, t)u s {- k - p - q, t) > . 

The tensors S , T, U are called cumulants of the second, third and fourth order, respectively. 
Following Orszag 16 , we obtain equations for S a fj(k,t) and T a/ g 7 (k, p, t). We have 

+ 2vk 2 )S afj (kp ) = 

k, ttt)Q aprj (p, k)T Pprj (- k, p, /) + f(p, -k, klt)Q fjpa { p, -k)T„ p(J (k, p, /)) (A.l) 

p 

and 

+ K fc2 + P 2 + g 2 ))r„, 3 7 (k, p, t) = -- ^(<5„ p<J (r, k)/(r, k, klt)Uf 3lp(T { p, q, r, t) + 

r 

Ql3p< r(r, p)f(r,p, klt)U aipa { k, q, r, /) + 

Q^fpa ( r ? q)/(D ?? Qt)Uaf3p<j{ k, p, r, t)). (A. 2) 

We then use a quasi-normal assumption and express U a fj pa ( k, p, r, t) in terms of S It implies 

( + K fc2 + P 2 + g 2 ))r„/3 7 (k, p, /) = 

2 

apa( P? k)/( Pj fc, (q? 0 H - Qo'p(j( q? k)/( ^ "f" 

Q (3 pa ( k, p)/( p, f^^o^k, ^)5*-y'(r(Q? Q (3 pa ( q? p)/(-q ? p> ^)^(k, 0 + 

^7 l od( k, (j) f ( k , 0^/5*Q-p(k, ^)^*/3 (t(p? H - Q 7/9<t( P? q)/(-P, ?, ^f)^a(j(k, /)^ P (p, *))• 

Before we proceed further, the Eddy Damping approximation must be introduced. For the 
non-rotating turbulence, it was found that the quasi-normal approximation alone is not sat- 
isfactory. In particular, the energy spectrum obtained in evolution would become negative 
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valued. Orszag 16 found that this unphysical behavior may be understood on the basis of 
improper relaxation time. Formal solution of T from its evolution equation introduces a 
memory integral. The memory integral should not involve only a viscous cut-off when the 
Reynolds number is large. Instead, the memory should cut-off because phase correlation 
between Fourier modes (or eddies) do not persist in the random convention fields of other 
eddies. The Eddy-Damped Quasi-normal approximation is proposed [Orszag 16 ]. The es- 
sential of the ‘Eddy-Damping’ is to approximate the fourth-order cumulants neglected by 
a linear damping term and to enhance the viscous term by p(k 2 + p 2 + q 2 ) + Pkpq- The 
parameter fj,k pq , which has the dimension of the inverse time, is the characteristic “decor- 
relation rate” of the third order moments. Lesieur 17 stressed that its choice is essential if 
one want to use the theory for quantitative predictions. For isotropic turbulence, Orszag 16 
recommended that 


f-^kpq — f^k T l-^p T f-Lq 


(A.3) 


where 


Hk ~ [k i E(k)] 1 l 2 


(A.4) 


is the inverse of the local eddy turn over time. This approximation is known as the Eddy- 
Damped Quasi-Normal (EDQN) approximation. 

As discussed in detail by Orszag 16 and Lesieur 17 , the EDQN approximation does not 
guarantee the realizability (the positiveness of the energy spectrum). A simple modifica- 
tion, called ‘Markovianization’, is needed. The essential of this approximation is that the 
exponential term varies at the nonlinear time scale r n i = [ v(k 2 + p 2 + q 2 ) + fJ,k pq ]~ 1 , which 
is much shorter than the time scale governing the evolution of ^ < uu >< uu > (on 
the order of the large-eddy turn-over time). The Eddy-Damped Quasi-Normal Markovian 
(EDQNM) approximation, therefore, assumes the effect of the time integration only over 
the exponential part memoral integrals: 


0 ! = 


exp {-[v{k 2 + p 2 + q 2 ) 


+ Vk P q](t - t )}/(— p, k, Dr)dr, 
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0 2 = f exp {-[v(k 2 + p 2 + q 2 ) + fi kpq \(t - T)}f(-q,k,ttT)dT 

Jo 

0 3 = [ exp {-[v(k 2 + p 2 + q 2 ) + Ii kpq \(t - T)}f(-k,p,ttT)dT, 

Jo 

6»4 = / exp {-[v(k 2 + p 2 + q 2 ) + /j, kpq \(t - T)}f(-q,p,ttT)dT 

Jo 

0 5 = exp {-[v(k 2 + p 2 + q 2 ) + Ii kpq ](t - T)}f(-k,q,ttT)dT, 

Jo 

0 6 = exp {-[is(k 2 + p 2 + q 2 ) + /i, kpq \(t - T)}f(-p,q,ttT)dT. 

Jo 


(A.5) 


Then we obtain 


Ta(3 7 (k,p,/) = i(0iQ apd{ P; k)N, 3 p (p, t)S 3 (J (q, t) + 0 2 Q apd{ *T k)5 / 3 ( 7 (p,i)5 , 7 P (q,i) + 

QoQfjpai-K p)S ap (k, t)S~ /a (q, t) + 04 Q [ 3 pa( — Q, p)5„ (J (k, t)S~ /p (q, t) +(A. 6 ) 
7 pa( k, q) N„ p (k, pS ^ p 0 qQ 7 pa( P? q)S aa (k,pSp p (p,p). 

Using (A. 6) we get expressions for Tfj pa ( — k,p,t) and T apcr (k,p,t ) in (A.l). Now we 
substitute these expressions in (A.l) to obtain 

(^ + 2 vk 2 )S afj (kp)= (A. 7) 

— 'y P f(Pi k) oLpaijPi k)(9iQ p^ rj ( — p, — k)Sp^(p^ t)S arj ((ii t) + 

v 

@ 2 Q {3£,t)( tyi k)S prj(Pi t ) S fj ^ ) T 
^3 Q p£,rj{k i p) S k,t)S ar] (q,t) T O 4 .Q QiP)$j 3 pi k^ t)S a ^Q^ T 

a £rj{k 1 j 3 ^( k , t)S pp (p, t ) T OqQ ( j£ ^ ( p, k , (p, ^)) 

— ^ P f(Pi —k-, Qt)Q p p a (p, — p, ^)5*p^(p, 0 T 

V 

@2$ a£rj( tyi k)S prj(Pi T 

@3Q p£p{. k*, p)Sa£,(ki T 64 . Q p£r}( Qi p) ^ ap{k 1 t) S T 

k , ^)-5* T ^6Q(7^?)( Pi Q) ^ ap{k 5 ^)-5* p^(p ? ^)) • 

The EDQNM evolution of the second order spectral tensor gives a complete description 
of the energy transfer and anisotropic structure of rotating turbulence. 

Taking trace in (A. 7) and combining terms, we obtain an equation for U(k,i). 

+ 2 vk 2 )U{k,t) = ^(a'(A:,p,g,fi,/)C(p,/)C(g,/)+ (A. 8) 
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b'(k,p , q , 0, t)U(k, t)U(q , /) + d(k,p, q, H, t)f7(£;, t)f7(p, /)). 


In (4.2) it is understood that k = p + q. The geometrically determined coefficients a'(k,p , g, 11, t), 
b'(k,p,q,Cl,t) and d(k,p,q,Cl,t) are given by 

a '(k,p,q,tt,t) = ~^f(p, k,klt)Q apa (p, k^Q^-p, -k)P pi (p)P ar] (q) + 

@2Q a(r]( — ~k)Ppr]{p)P(j({ ( i)) ~ 

—k? klt)Q a pa(Pi ~k)(0\Q a^p^—p, k)P p ^(p)P crr] (q) + 

02Q a ^ n {-q,k)P p , n {p)P ai {q)) (A. 9) 


1 

4 


b (k,p, q, H,t) — f {p^k ,klt^Q ap a{p^k)(9 3 Q p^pik , p^)P a ^{ k)T > (J ^(q) + 

@aQ p^r]( — P)Par]( — k)T > (J ^(q)) — 

f(p, -k,ttt)Q apa (p, -k)(d 4 Q piv (-q,p)P av (k)P ai (q) + 0 3 Q piv (-k,p)P ai (k)P av (qj^.W) 


c (k,p, q, H,t) — — f(p,k, klt^Q aprj (p, k)(9 3 Q a ^ p (k^ q)Pa^{ ^P)Ppp (p) + 

@&Q <j(r](— Pi <l)Pari( — k)Tp^(p)) — 

— —k? klt)Q ap a(Pi ~k){0 3 Q^^, p ( y —k^ q)P a ^(k)P pr] (p) + 

@&Q <j(r](—Pi < l)Par]0^)Pp({ P))- 

We note that the geometrically determined coefficients a', d and d contain / and 9. 
The coefficients 9 are given by (A. 5). 
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